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Figure Captions 

Figure 1 The ratios (2.6) of the dynamical structure factors at p = tt versus the 
logarithm of the euclidean time r: a) longitudinal 7 = 7r/2; b) longitudinal 
7 = 0.37T. The dashed curves represent the estimate of the thermodynamical limit 
according to (2.10). The dotted curves are the integrals (2.13) 

Figure 2 Same as figure 1 for the isotropic case 7 = 0. 

Figure 3 Same as figure 1 for a) transverse 7 = l/37r; b) transverse 7 = 1/2tv. 
Figure 4 The finite size behaviour of lnr along the horizontal lines Ri = 1/4, 1/2, 3/4 
shown in figure 3(b) 

Figure 5 The finite size behaviour of ln(r + t ),t = 0.2 along the horizontal lines 

R = 1/4, 1/2, 3/4 shown in figure 2. 
Figure 6 The ratios (4.1) of the dynamcal structure factors at noncritical momenta 

versus lnr: a) isotropic p/ir = 1/4, 3/4, N = 8,16; b) isotropic p/ir = 1/2, N = 

4,8,12,16 

Figure 7 The quantity p(t,p = vr/2,7 = 0,iV) - defined in (4.3) - versus t - defined in 

(4.4). The dashed curve is the prediction (4.5,6) 
Figure 8 The ratio (4.1) of the dynamical structure factors at fixed euclidean time 

r = 0.5 versus momentum p and N = 4, 6, 16 
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Appendix A. Excitation Energies and Transition Probabilities for N = 16 

In the following table we present the excitation energies u n = E n — E and the 
transition probabilities t n = | (n|s 3 (0) |0) | 2 for the isotropic case at fixed momenta 
P /tt = 1/4,1/2,3/4,1: 
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can be described adequately by the ansatz (2.10). The resulting thermodynamical 
limit (2.8) for the ratio (2.6) is in good agreement with the exact result of Katsura 
et al.[4] at 7 = 7r/2 but deviates from the prediction (2.13) of Miiller et al. [10] for 
smaller values of 7. 

(iii) At the critical momentum p = tt finite size effects are large in the transverse case 
and increase rapidly with increasing values of 7. It was demonstrated in section 3 
that this behaviour signals the emergence of the (nonintegrable) infrared singularity 
in the tranverse structure factor. 

Therefore we can conclude that the euclidean time representation is particularly 
suited for the study of finite size effects, which allows for a crude estimate of the 
thermodynamical limit. This estimate is already precise enough to check the gross 
features of a model ansatz [10] for the dynamical correlation functions in the spectral 
(cu)-representation. To resolve the fine structure, however, - i.e. the detailed form and 
cutoff of the nonsingular contributions- one has to know the euclidean time dependence 
very precisely. In other words: The reconstruction of the spectral (cu)-representation 
from the euclidean time (r)-representation demands for an anlytic continuation and 
tiny errors in the r- might produce large errors in the cu-representation. Of course this 
problem would not occur if finite size effects could be analyzed systematically in the 
spectral representation. The continued fraction method used in Reference [14] might 
open this possibility. 

Our exact results on small systems might be useful as a test for those, who plan 
to develop approximative calculations on larger systems. For this purpose we present 
in appendix A our results for the first nine excitation energies with the corresponding 
transition probabilities at fixed p/ir = 1/4, 1/2, 3/4, 1; N = 16; 7 = 0. 
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where 

tj = t 1 -^ exp^r). (4.4) 

According to (4.2), pj(tj,~f,N = oo) should be linear in tj for tj — > oo. In figure 7 we 
show this quantity for the isotropic case at p = tt/2 and N = 4,8,12,16. The linear 
behaviour in t becomes more and more apparent with increasing system size. We can 
compare our results for the longitudinal case with the prediction: 

Mr,P,T = 0,7) = J 3(T, !! ,7) v (4.5) 
I 3 (t = 0,p,7) 

h(r,p,i) = f M" 2 -"l)- a3 (" 2 2-" 2 )- 1/2+a3 e-" T (4-6) 



which follows from the ansatz of Miiller et al.[10]. In (4.6) the lower and upper 
integration bounds are given by (1.3) and: 

2tt p 

U2{P,V = — sin 7 sin-, (4.7) 

7 z 

respectively. The prediction is found between the finite system results for N = 12 
and N = 16 (cf. the dotted curve in figure 7). Better agreement with the expected 
behaviour in the thermodynamical limit can be achieved for example by lowering the 
upper integration bound cu 2 in equation (4.6). This, however, would violate the sumrules 
of Ref. [13]. To our knowledge there does not yet exist a generalization of the ansatz of 
Miiller et al. [10], which respects the sumrules in Ref. [13]. 

In figure 8 we show the ratio R as fuction of the momentum p (0 < p < it) and at 
fixed: r = 0.5, ~{/tt = 0.0 N = 4, 6,..., 16. All these data points nicely follow one 
curve in the 'scaling'-variable p. Again this means that the thermodynamical limit is 
seen already on small systems. 



5. Conclusions 



In this paper we started a first attempt to extract the dynamical structure factors of 
the XXZ-model from a complete diagonalization of the Hamiltonian on finite rings with 
N = 4,6, ...16 sites. We studied the normalized ratios (2.6) of the dynamical structure 
factors as function of the euclidean (imaginary) time r and found the following features: 

(i) Away from the critical momentum p = tt finite size effects are small except for the 
large r-limit, where we find a clean signal for the threshold singularity (1.4). 

(ii) At the critical momentum p = tt finite size effects are still small in the longitudinal 
case at 7 = 7r/2 but increase for decreasing values of 7. These finite size effects 
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So far our discussion of the nonintegrable infrared singularity is restricted to the 
transverse case with anisotropy 7r/2 > 7 > 0. In the isotropic case 7 = 0, a\ = a 3 = 0.5 
the integral (3.3) diverges logarihmically for x — > 0: 

I^x,^ = 0) = -i(oo, 7 = 0)lnx (3.8) 

From (3.8) we predict that the 'half width' in the isotropic case: 

b[r(i?i,7 = 0,7V) + T ] = (i?x - l^nu^N) + C{R 1 ) + ... (3.9) 

increases with IniV but with a slope (1 — i?i) depending on R\. Again this behaviour 
is visible even on small systems, if we make a proper choice for r = 0.2, as can be seen 
from figure 5. The observed slopes have the correct i?i-dependence. Their absolute 
value differs from the expectation by about 10%. 

4. Dynamical Structure Factors in the Noncritical Regime 

Leaving the critical momentum p = tt the threshold singularities in the dynamical 
structure factors change in position and strength. The singularity moves according to 
(1.3) to nonvanishing frequencies u = iO\(p } 7); its strength is reduced to 
{u — LOi)~ aj . The threshold singularity is integrable now, since ctj < 1 for j = 1,3 and 
the Laplace transforms of the transverse and longitudinal structure factors exist for all 
nonnegative r- values, if < p < tt. The ratios: 

rm-p-~-n) = ,f ,T -f;°-i N L ("> 

Sj(T = 0,p,l = 0,7,iV) 
converge to a limiting function in r. Its large r-behaviour is given by the threshold 
frequency (1.3) and the strength of the threshold singularity: 

R 3 (r,pn) ™ exp(-cu 1 T)r^- 1 . (4.2) 

We have determined the ratios Ri and R 3 as function of r at fixed noncritical momenta 
and for various values of 7. We found extremely small finite size effects. Moreover 
the ratios (4.1) almost coincide for the different values of 7 and for the longitudinal 
and transverse case. We only observe a weak dependence on the momentum p. As an 
example we present in figure 6 the results for the isotropic case (7 = 0) at fixed momenta 
p/tt = 1/4,3/4, 1/2. The first two momenta can be realized for N = 8, 16. Here finite 
size effects cannot be resolved in the plot of the ratio (4.1). The origine of the lnr-axis 
has been shifted by 2 in order to present the results for p = 7r/2, N = 4, 8, 12, 16. Again 
the ratios (4.1) coincide for N = 8, 12, 16, whereas the result for N = 4 is found left to 
them. 

Significant finite size effects only appear for large values of r and are visible in the 
quantities pj(tj,~f, N) which are related to the ratios (4.1) via: 

R j (r, 7 ,N) = [l + p j (t j , 7 ,N)]-\ (4.3) 
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The second one is assumed to be free of such a singularity. We have introduced an 
exponential cutoff for the high frequency contributions in the first term with a parameter 
r , which will be fixed below. The factor A(u/u>i, 7) is supposed to describe the 
approach to the infrared singularity. For u\ — > 0, A(lo/uji , 7) is assumed to converge to 
a nonvanishing value for the residue ^(00,7) of the infrared singularity. Starting from 
(3.1) we find for the ratio (2.6) in the limit u\ — > 0: 

^(M,^ 00) = /l( " l(T + To) - 7) + "\" 1/{(T,7) , (3-2) 

17 } / 1 (c 1 ro,7)+^ 1 - 1 /{(0, 7 ) 1 } 

where: 

00 

Ji(a, 7 ) = j dytf-iy^e-^Aiy^) (3.3) 
1 

and 

00 

/((r, 7 ) = J dioS'^^e—. (3.4) 


In the combined limit: 

N — > 00, r — > 00, x = u\(t + r ) fixed, (3-5) 

we expect the ratio (3.2) to converge to a scaling function I\(x } 7)// 1 (0, 7), since 
2«i — 1 > for < 7 < 7r/2. The scaling function (3.3) depends on the parametrization 
of the first contribution in (3.1). The cutoff parameter r enters in the finite size 
corrections to the scaling variable x. The scaling curve depends explicitly on A(y,~f). 
The scaling behaviour of the ratio (3.2) in the combined limit (3.5) has an immediate 
consequence for the 'half width' : 

ln(r(i? 1 , 7 ,iV) + To) = - lni^iV) + In x(i? l5 7) + ... (3.6) 

which diverges for N — > 00. Therefore on finite systems, the signature for the emergence 
of the nonintegrable infrared singularity is a linear increase of the half width with In N 
and with slope 1. This behaviour should be observable not only for Ri = 0.5 but for 
any fixed value of Ri between and 1. In figure 4 we have plotted the left hand side of 
(3.6) with r = versus — Inuji(N) for i?i = 0.75,0.5,0.25 and 7/^ = 0.5. The linear 
behaviour in —InuJi(N) is clearly seen and the slopes 0.82,0.92,0.96 are found to be 
rather close to the expectation, namely 1. The second term in (3.6) (x(i?i,7)) is the 
inverse of the scaling function I\(x } 7)// 1 (0, 7). The behaviour of the scaling function 
for small values of the scaling variable x = u>\(t + t ): 

h(x -+ 0,7) - /!(0,7) = -i(oo, 7 )r(l - 2a 1 )x 2a ^\ (3.7) 

is governed by the exponent «i of the infrared singularity. 
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curves in figures. 1(a), (b) represent the resulting scaling curve i? 3 (r, 7) for 7 = 0.57T 
and 7 = 0.37T, respectively. The exponent /3 3 (i? 3 ,7) decreases with decreasing values 
of R 3 and 7, which is a signal for increasing finite size effects. This feature is easy 
to understand. With increasing values of r the low energy excitations get a stronger 
weight in the Laplace transform (2.2). On the other hand, the spectrum of low energy 
excitations is particularly sensitive to the finiteness of the system. We can compare our 
scaling curve i? 3 (r, 7) with the ansatz of Miiller et. al. [10] for the longitudinal structure 
factors: 

2A 



S 3 (us,p = 7r,T = 0,7,iV = 00) 



B(l -a 3 ,l/2 + a 3 ) 

0(^2(7) - u) 



CU 2 «3(cu 2 ( 7 )2 -^2)1/2 



[2.11] 



where B(x, y) is the betafunction. In this ansatz the high frequency excitations are cut 
off at 

, , sin 7 

co 2 ( 7 ) = 2tt '-. (2.12) 

7 

The ansatz coincides with the exact result of Katsura et al.[4] in the XX-limit 7 = tt/2 
and leads to the following expression for the ratio (2.6): 

1 

Mr,l) = W n? 2 1/9 ^ , [dxx-^(l-x*r-^e-^ (2.13) 

#(1/2 - a 3 , 1/2 + a 3 ) J 



The integral (2.13) is represented in figures 1(a), (b) by the dotted curves. For the XX- 
case 7 = 7r/2, we find good agreement of our determination of i? 3 (r, 7) with the exact 
result (2.13). This agreement supports our hypothesis, that the finite size effects can 
be parametrized adequately by the ansatz (2.10). For 7 = 0.37T, our finite size analysis 
leads to a scaling curve i? 3 (r, 7) which differs significantly from the prediction (2.13). 
This discrepancy might originate from the sharp cutoff (2.12) in the high frequency 
excitations. 



3. The Nonintegrable Infrared Singularity in the Transverse Structure 
Factors at p = it and T = 

Let us assume that the transverse structure factor can be split into two parts. The first 
one contains the nonintegrable infrared singularity ((1.4) for j = 1): 

S^p = tt,T = 0,7, W ^ 00) = e( f_~Tl e—°A(ulu ir t) + 

[CO L0 1 ) 

Si(u,-y). (3.1) 



between the groundstate energy E and the energy E\ of the first excited state. 

In the following we consider the Laplace transform (2.2) normalized to the corresponding 

static structure factor: 

z? ( at\ S j( t ,P = k,T = 0,7,7V) 

Rj{T,l,N = 7^7 J: =, J: ttt, J = 1,3. 2.6 

^•(r = 0,p = 7T,T = 0,7,7V) 

By construction these ratios are monotonically decreasing with r and varying between 
1 and for r > 0. Figures l(a),(b);2;3(a),(b) present the ratio (2.6) for the longitudinal 
case at 7/7T = 0.5, 0.3, the isotropic case 7 = 0, and the transverse case at ^/tt = 0.3, 0.5, 
respectively. Going from figure 1(a) to figure 3(b) we observe the following characteristic 
features: 

(i) Finite size effects are small in the longitudinal case at ^/tt = 0.5 and increase with 
decreasing values of 7. 

(ii) Finite size effects are large in the transverse case and increase with increasing values 
of 7. 

(iii) The 'half width' t(R 1 = 0.5,7)- defined as the value of r, where the ratio (2.6) has 
dropped to Rj = 0.5- moves systematically: 

r(i? 3 = 0.5,7 = 0.5tt) < t(R 3 = 0.5, 7 = 0.3vr) < t(R = 0.5, 7 = 0) 
< T (R 1 = 0.5,7 = 0.3tt) < t(R 1 = 0.5,7 = 0.5tt). (2.7) 

As will be shown below, this property is related to the strengthening of the infrared 
singularity (1.4) according to (1.5): 

a 3 (0.57r) < a 3 (0.37r) < a(0) < ai(0.37r) < ai(0.57r). 

The integrability of the infrared singularity in the Laplace transform (2.2) for the 
longitudinal case means that the ratio (2.6): 

i? 3 (r, 7 ,iV^^>) = i? 3 (r, 7 ) (2.8) 

converges to a scaling curve i?3(r, 7). The large r-behaviour of the scaling curve is given 
by the infrared singularity (1.4): 

^(m)™^" 1 . (2.9) 

We have tried to determine the inverse r = r(i? 3 ,7) of the scaling curve from a finite 
size analysis of our results for N = 4, 6, 16 with a parametrization: 

lnr(i? 3 ,7,iV) = lnr(i? 3 ,7) + (2-10) 

We fixed the parameters r(i? 3 , 7), B(R 3} 7), /3(i?3, 7) at N = 12, 14, 16. We checked the 
validity of (2.10) by comparison with our low N data (N = 4,6,8,10). The dashed 
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The outline of the paper is as follows. In Section 2 we present our results on the 
dynamical structure factors at p = tt and N = 4, 6, 16. We make use of the euclidean 
time representation, which is particularly suited to analyze finite size effects. In Section 
3 we demonstrate how infrared singularities emerge on finite systems. In Section 4 we 
report on our results for noncritical momenta p < tt. 



2. Finite Size Behavior of the Dynamical Structure Factors at p = it and 

T = 



In the critical regime: 



p 



7T, 



T -> 0, N 



oo, 



(2.11 



the dynamical structure factors (1.4) develop infrared singularities. We want to 
investigate how these singularities show up on finite systems, where the structure factors 
are sums of ^-function contributions. Finite size effects are not so easy to analyze and we 
therefore look for a smoothening procedure which allows to extract the thermodynamical 
limit from finite systems. For this purpose let us consider the Laplace transforms of the 
structure factors. At T = } p = tt they acquire the following form: 



Sj(r,p 



TT 



r = o, 7 ,iv) 



n>0 



Po 



tt e 



-t(E„-E )\ 



W*;(0)|0)|' 



(2.2) 



The variable r can be interpreted as an euclidean time. At r = we recover the static 
structure factors which behave for N — > oo as [12]: 



S 3 (t = 0,p = 7r,T = 0,7,7V 



oo = r,; 



(7): 



n ">fc(7)-l 



N 



i-v 3 (i) 



with critical exponents [9]: 

^1(7) = ? ?3(7)" 1 



1-1. 

7T 



(2.3) 



(2.4) 



According to (2.3,4) the longitudinal structure factor stays finite whereas the transverse 
one diverges. This divergence originates from the infrared singularity (1.4) which is 
integrable for the longitudinal case (20:3(7) < 1) but nonintegrable in the transverse 
case (20-1(7) > !)• On finite systems the infrared singularities are not visible directly 
due to the gap: 



u^p = 7T, 7 ,iV) = E 1 -E = OiN- 1 ) 



(2.5) 
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1. Introduction 

Based on a complete diagonalization of the XXZ-Hamiltonian: 

N 

H = 2^2[s 1 (x)s 1 (x + 1) + s 2 (x)s 2 (x + 1) + cosjs 3 (x)s 3 (x + 1)] (1.1) 

x=l 

on finite rings with N = 4, 6,..., 16 sites and anisotropy parameter ^/tt = 

0. 0, 0.1, 0.3, 0.4, 0.5 we reported in References [1] and [2] on the thermodynamics and the 
static structure factors at finite temperature. In this paper we continue our numerical 
investigation of the XXZ- model with an analysis of the dynamical structure factors: 

Sj(u,p,T,f,N) = Z~ x - (E n - E m ))8(p- p n + p m ) x 

KnlsMHW J = 1,3. (1.2) 

Z is the partition function and |n) denotes an eigenstate of the Hamilton operator and 
of the momentum operator with eigenvalues E n and p n , respectively. The dynamical 
structure factors contain the information on the transition probabilities | (ra|sj(0) \m) | 2 
between the eigenstates n and m with an excitation energy u = E n — E m and a 
momentum transfer p = p n — p m . At T = and p = tt the XXZ-model is known to be 
critical and here one expects that quantum effects become most important. Therefore, 
most of the previous studies were concentrated on the grounstate behaviour. There exist 
analytical results on the dynamical correlation functions for the special case 7 = 7r/2, 

1. e for the XX-model [3-6]. This model can be mapped on a free fermion system [7] 
by means of a Jordan Wigner transformation. For the general case < 7 < 7r/2, the 
spectrum of the lowlying excitation energies has been exploited by des Cloiseaux and 
Pearson [8]. In particular it was found from the Bethe ansatz solution that there is a 
lower bound 

7T 

E n - E m > ^1(^1,7) - E (p ,7) = ^i(p,7) = ~ sin7sinp (1.3) 

7 

for the excitation energies depending on the momentum transfer p = pi — p . 
Approaching the boundary (1.3) the structure factors diverge [9]: 

5>,p,T = 0, 7 ,iV = 00) = A { ^~^l y (1.4) 

where 

«i(7) = ^(l + 2 ), a 3 ( n ) = ^^. (1.5) 

Z V 7T/ 7T — 7 

Starting from these rigorous results G. Miiller and collaborators [10] made an ansatz for 
the dynamical structure factors which has been applied successfully on the description 
of neutron scattering data [11]. 
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